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o3 ■ Abstract 



Using a probabilistic interpretation of resonant states, we propose a formula useful to calcu- 
late the lifetime of a resonance using square-integrable real basis-set expansion techniques. Our 
approach does not require an estimation of the density of states. The method is illustrated with 
calculations of s and p resonant-state energies and lifetimes. 
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I. INTRODUCTION 



Resonant states are of great importance in quantum physics since the pioneering work 



of Gamow modeling the a decay in radiative nuclei [1[. Resonant states appear in nuclear 



QE 



nn 



physics [l|, |2|, atomic and molecular physics [3|, |4[, and more recently, in nanophysics 



a 



01?]. 



8MlO|| presented a probabilistic interpre- 



In a series of papers Hatano and collaborators 
tation of resonant states. In these works, a resonance is defined as an eigenfunction of a 
Hamiltonian with Siegert boundary conditions, that is, when the potential goes to zero, the 
wave functions is only an outgoing wave ll|. Due to Siegert conditions, the problem is 
not Hermitian; the eigenvalues are complex, the wave functions are not square integrable 
and the particle number around a central volume exponentially decays in time because of 
momentum leaks from this volume. 

Many methods were developed in order to calculate the complex eigenvalues of resonant 
states. The most widely used methods are based in complex scaling transformations that 



turn the resonance state into a square- integrable function 12J, ll3[ . In particular, a modifica- 
tion called exterior complex scaling |l4[ was used recently in many photoionization problems 



15 



16|. 



However, due to its simplicity, there also exist several techniques to calculate the real 
part of a resonance eigenvalue using the Ritz-variational method with a square integrable 
real basis-set. These methods, called in general stabilization methods 17|], use the abrupt 
change of a physical quantity when a variable is "crossing" a resonance. In particular 



different authors used the variational energies 18 



properties as the double-orthogonality condition 20] to obtain the real part of the energy. 



19| . the von Neumann entropy [2J, or 



0, 



In order to calculate the imaginary part of the resonant energy, all these methods use an 
app i 

3, 



approximation of the density of states fitting the numerical data using a Lorentzian function 
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211 • 



In this paper we use the formalism of Hatano |8NlO| to obtain an expression of the 
imaginary part of the energy useful to calculate it using only a stabilization method, without 
an extra fitting or approximation. 

The paper is organized as follows. In Section II we use the Hatano formalism to obtain a 
relation between the real and imaginary part of the energy and the density of probability in 
the central region. In Section III we apply the results of the preceding section together with 



real basis-set expansions to calculate the complex energy of a resonance. Finally, Section IV 
contains the conclusions with a discussion of the most relevant points of our findings. 

II. DEFINITION AND PROPERTIES OF RESONANCES 

Here and elsewhere we use atomic units, H — 1 ;m — 1. A general solution of the time- 
dependent Schrodinger equation 



d 
Hip(x,t) = i—ip(x,t) 



(1) 



obeys a continuity equation 



dp(x,t) 
dt 



+ V • J(x, t) = 



(2) 



where p and J are the standard density and current of the probability 22] 



p(x,t) = 
Eq. ([2]) could be written in integral form 



x,t)\ 2 ; J(x,t) = -[iP*(x,t)V^(x,t)-^(x,t)Vr(x,t)} . (3) 







p(x, t) d d x 



J(x, t) ■ dS . 



(4) 



9t jq jqq 

where d is the spatial dimension, Q an arbitrary volume, and dQ its frontier. Following 



Hatano 



-|l0|]. we define Nq, the number of particles inside a volume Q 



N u (t) 



From the definition of J 



p(x, t) d d x . 



(5) 



22|, Eq. (EJ) takes the form 



J^tfn(t) = ~ Re [l r&t)pip(x,t)-dS) . 



(6) 



This equation, that express the particle-number conservation inside the volume Q, corre- 
sponds to Eq. (15) of reference [9(. If the wave function goes to zero faster than r"^ 1 ) 
for large values of r the RHS of Eq. ([6]) goes to zero, expressing the conservation of the 
normalization, as happens for bound states. Resonant eigenstates are not Hermitian solu- 
tions of the Schrodinger equation, and the wave functions have exponential divergences, Eq. 



([6]) describes a flux of particles outside any volume Q, even in the limit Q — > Mr, and the 
particle-number is not conserved. 

Hatano suggested how to maintain the probabilistic interpretation of the wave function 
for resonant states |8l-ll0|. The resonance is interpreted as a metastable state, which is 
localized inside a volume Q at t = 0. Hatano defined a time dependent volume Q(t) by the 
condition 



j t N m (t) 



0. 



(7) 



The reasonable initial condition for this equation is iVj2(o)(0) = 1. This condition together 
with Eq. ([7|) imply that N^t)(t) = 1 Vt > 0. The expectation value of an operator O is 
calculated inside the volume Q(t) as 



(6) 



U(t) 



(1> 6 4>)cm 



(^ I V»)n(t) 
which is well defined for all times. 

Eqs. ©,06]), and (|7|) give an equation for dQ(t) 



(8) 



Re 



( dx 




0. 



(9) 



In order to obtain solutions of Eq. (jUJ), we have to particularize the system. We restricted 
our study to solutions of the time- dependent Schrodinger equation with fixed energy 



ip(x,t) 



-iEt 



iPe(x) 



(10) 



where E and i/je{x) are the eigenvalue and eigenfunction of the time- independent Schrodinger 
equation. We assume one-particle Hamiltonians with central potentials that tend to zero at 
infinity. Then, for central potentials we can use the reduced radial Schrodinger equation for 

/-waves 



H t u Ej i(r) = Eu E ,i(r) 



(11) 



where 



Hi 



+ ■\rrr Z + V(r) , and ijj e {x) 



u E,l( r 



2 dr 2 



2r 2 



Y hm {n) 



(12) 



where Y/ >m (f2) are the spherical harmonics 22) . 



Resonant states obey non-Hermitian Siegert boundary conditions |ll| . u E) i(r) ~ e lkr for 
r — )■ oo, where fc = y2E. Siegert states are non normalizable, so they do not belong to the 
Hilbert space, and their complex eigenvalues E = £ — zT/2 are interpreted as energies E 
and inverse lifetimes V of metastable resonant states. 

By symmetry, the solution of Eq. ()7j) for the volume Q(t) is a sphere of radius R(t), 
Q(t) = B(R(t)), with the initial condition Rq = R(t = 0). In this case, Eq. ([9]) gives the 
evolution of R(t) 



R(t) = Im l ' drUEAr] 



(13) 

r=R(t) 



U E,l( r ) 

The RHS of this equation does not depend explicitly on t, and the formal solution is 

t = / t r- . (14) 

\ UE,i(r) J 

Even the resonant-states functions are not square-integrable, the initial condition for Ng 
fix the arbitrary normalization constant 

/[•Rq 
\iP(x,0)\ 2 d 3 x = I \u Ejl (r)\ 2 dr = 1. (15) 

With this condition, the expression for NB(R(t)) (t) takes the form 

r R{t) 
N B (R(t))(t) = e" rt / dr\u E ,i(r)\ 2 = 1, (16) 



and the derivative with respect to t gives 



f* t] dr\u E ,i{r 
combining this equation with Eq. ffT3|) we obtain 



«'W»'' R(t), (17) 



Im [u* E>l (R(t)) d r u E ,i{r)\ r=m ^ 

J dr\u E ,i{r)\ 
This simple expression for T is not convenient for numerical calculation with real basis-sets. 
With this purpose in view, we write the explicit expression for the asymptotic behavior of 
the wave function 



u E ,i(r) = Ce tkr v L (k,r) 



(19) 



where C is a normalization constant given by Eq. ( TT5l) . Eq. (TT91) in Eq. ffTBl) gives 



Im < ik 



d r vi(k,r) 



Vi(k,r) 
By definition k 2 = IE = 2£ — «r, which gives 

Im(i k) = Re(k) = 



\u E ,l(R(t))\ 2 

-mi J J* (1) dr\u E ,i(r)f 



<£2 



£ 



1/2 



Finally, using Eq. ( 12~T|) in Eq. (I2"0"j) . we arrive at the expression 



(20) 



(21) 



£ + \ £ 2 



+ Im 



vi(k,r) 



r=R(t) 



\uE,l(R(t))f 



L i dr\u E ,i{r) 



For the important case treated in references 



(22) 



9] of potentials with a finite support, 



V(r) 



V(r) if r < To 
ifr > r 

we have an explicit expression of vi for Z-waves, valid for r > r 



(23) 



v t (k,r) = 5^( _1 



3=0 



,i V + jV- i 

j!(Z — j)! (2ifcr)* 



(24) 



In particular, for s-waves, I = 0, Vo(k,r) = 1, and then Eq. (1221) reduces to a linear relation 
for r 2 , which gives 



2£ + 



\u E ,o(R(t))t 



1/2 



Ko(fl(t))| S 

-fl(t) 



J dr|u Sj0 (r)p 



(25) 



.2/* (t) dr| M£ ,o(r) 12 

Equation (I2"2"j) . and for / = 0, Eq. ( I2"5"j) relate T with real magnitudes, £ and \ue,i(R)\ • This 
fact makes these equations useful tools to calculate T using square integrable real basis sets, 
as we show in the next section. 
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III. THE RITZ- VARIATIONAL METHOD AND RESONANCES: NUMERICAL 
EXPANSIONS 

In this section we use Eqs. (j22p and (I25p to calculate the imaginary part of the energy 
of a resonant state, applying the Ritz-variational method with real square-integrable basis 
sets. 

We exploit three facts of the variational expansion: 

(z) There are several accurate methods to calculate the real part of the eigenvalue, or 
energy, £, of a resonant state. 

(ii) The variational method gives good approximations to the exact (non-normalizable) 
densities p(r) where the resonant states are localized (see figure [2]). 

(Hi) Eq. (1221) involves just real quantities that could be evaluated with Ritz-variational 
wave functions. 

Even Eqs. (1221) and (1251) are valid for R > r , because item (ii), in the numerical 
calculations we take R = R = r , and then J ° dr \v,E,i(r)\ = 1. 

A. 1 = 

We begin with the simple case of s— waves, I = 

For a clear notation, we will omit the subindexes E and 1 = 0, then u(r) for an arbitrary 
potential V(r) has the form 

(ii < (r) ifr < TV, 

. ■ (26) 

u < (r )e jfc(r " r °) ifr > r 

Then, in this case Eq. (|2^|) takes the form 



2g+ i^ii |w < (ro)|2 _ (27) 




Eq. ( 12T|) involves two real quantities, |M < (r )| and £, both quantities are well approximated 
by applying the Ritz method using a real square-integrable basis set truncated at order 
N, {§i}i ■ In this approximation, the Hamiltonian Hi is replaced by a N x iV Hermitian 
matrix [ifj]jj = (&i\Hi\&j), and we obtain N eigenvalues E n and eigenvectors a^ n \ The 
corresponding orthonormal eigenf unctions are 



2V 



€ N \r) = Y,^^ 



N 
V^ / (n) 

a; 



i=l 



2=1 



1 ; n 



JV. 



(28) 



In particular, we used the double-orthogonality method (DO) [20j in order to calculate the 
real part of the resonant energy 8^ and the variational square-integrable approximation to 
the resonant wave function vS N '(r). This method assumes that the potential depends on a 
parameter A, and when A is varied on an interval [Ax, \r\ a given eigenvalue n$ crosses the 
resonant energy value at A no , as illustrated in figure [TJ The method uses the fact that in 
both sides of the interval, A^ and \r, the n eigenvalues and eigenvectors in the left and in 
the right of the avoiding-crossing zone correspond to different states of the quasi-continuum 
and the resonant state is orthogonal to both of them. We define the double-orthogonality 
function 



D n {\) = |(^(A L ),^(A)>| 2 + |(^(A fl ),^„(A))| 2 , for \ L < A < \ R . 



(29) 



Because the eigenf unctions are normalized, < D n {\) < 2. For a given eigenvalue no, 
we define the localization of the resonance A no as the value of A where D no (X) reaches its 
minimum, that is, where the eigenfunction has a minimum projection onto to the quasi- 
continuum states, at it is shown in figures |3] and [6j This method has the advantage over 
other stabilization methods in that we have to solve the variational problem just a single 
time. The price we pay is that A n is also an output of the method and we cannot choose it 
arbitrarily. The best approximation of the resonance is defined as 



A 



— min D no (\) ; £(A no ) — E no (X no ] 

Xe[X L ,X R \ 



(30) 



Once we determine the optimal A„ , its eigenfunction iph (r) |20| together the normal- 
ization condition Eq. (IT51) give 



u iN) (r) 



0> 



no \ 



ELaf"* 



where 



Jo" ^oV) 



kAR) 



dr 



^N (TV) {N) 



KU <X v ^>o) 



R 



$m(r)$„(r)o?r. 



Then, from Eq. (|2"TJ) we obtain for r^ 



(31) 



(32) 



L ' I ' ^ 
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— i- 
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FIG. 1. (color online) (a) The first 30 eigenvalues of the N = 100 Hamiltonian matrix of the I = 
block of the potential Eq. ()34|) as a function of the barrier height A. (b) Same as (a) plus exact 
(green line) and approximate (blue dots) energies of the resonant state. 



r (JV) 



Vff (r ) 



N 



i,3 



(N)(N) 



=1 ">i 



a- a 



T'k^o) 



rltiP{r ) 



\ 



-,1/2 



V 



2£ 



N (N) (JV) 

tj=l a i a j 



kj( r o) 



+ 2£ 



(33) 



As a particular case, we calculated resonant states for an exactly solved problem, the 
well+barrier potential 23 ] 



V(r) 



-Vc, ifr < A 

-A if A < r < r 



(34) 



where all the parameters are positive. The exact wave functions are different combinations of 
exponential functions in each sector, with continuous logarithm derivative at r = A and r = 
tq. The exact energies for bound, virtual and resonant states are given as solutions of three 
different transcendental algebraic equations, which are obtained applying the corresponding 
boundary condition at r = r 23]. 

The calculations were done as a function of the barrier height A, with fixed values of 
Vq = 0.15, A = 5, and ro = 6. A convenient orthonormal basis set is given by 



P 0.2 - ■ 




0.005 - 




FIG. 2. (color online) p{r) for a I = 0-resonant state, exact (black line) and Ritz-DO approximation 
with N = 100 (red line) of the potential from Eq. (JMJ). (a) global view, (b) Region < r < r$ = 6, 
where the resonant state is localized, and (c) The exterior region r > r$ = 6, where the resonance 
diverges and the variational expansion is a stationary wave. 



D 



n 




0.25- 



FIG. 3. D n for the I = block of the potential Eq. (|34|) as a function of the barrier height A 
for N = 100 and n = 2, . . . , 30. The minimum of each curve is defined as the localization of the 
resonant state. 
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FIG. 4. (color online) Exact (red line), N = 100 (turquoise dots) and N = 500 (black squares) 
values of T against the barrier height A for the I = block of the potential Eq. ([34 



$;(r) 



y/(i + l)(i + 2) 



r/2 L?Hr) ; i 



1 N. 



(35) 



where L\ (r) is the Laguerre polynomial of degree i and order 2 [24]. 

In figure [Tfa) we show the first thirty eigenvalues given by the Ritz method with N = 100. 
In figure QJb) we add the exact resonant energy curve and the values calculated with the 
DO n functions, which are shown in figure [3] for N = 100, and n = 2, . . . , 30. Once the 
energy was obtained, the width T of the resonance is calculated using Eq. ( |33l) . In figure 
H] we show r(A n ) for two different sizes of the basis set: N = 100, for n = 2, . . . , 30, and 
N = 500, for n = 2, ... , 140. Our data show an excellent agreement with the exact curve 
r(A), that is also included in the figure. 



B. 1 = 1 



From Eq. (1241 . the function vi(k,r) for p- waves takes the form 
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0.04 = 




FIG. 5. (color online) The first 30 eigenvalues of the N = 100 Hamiltonian matrix of the I = 1 
block of the potential Eq. (|34p as a function of the barrier height A. The approximate energies of 
the resonant state calculated with DO are also shown (blue dots). 



vi(k,r ) =l + i — • 

kr 

In this case, it is convenient to re-write Eq. ( l22l) as a function of x = Im(k) 



(36) 



x 3 + 



j_ N r o)r 

r 2 



x 



£ 



K r o)r 



2rl 



x 



£ 



K r o)r 







(37) 



2r 

The definition of resonance establish that Im(k) > 0, and we proved for each case that we 
studied that Eq. (1371) has an unique positive root. Finally T is obtained from the definition 
of k as 



T = -2x p ^2S + x 2 p , (38) 

where x p is the positive root of Eq. (1371) . 

We calculated the inverse lifetime T for a / = 1 resonant state of the potential Eq. ( )34|) 
as a function of the barrier height A, with fixed values of V = 0.3, A = 5, and r = 6. 
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FIG. 6. D n for the I = 1 block of the potential Eq. (|34j) as a function of the barrier height A 
for N = 100 and n = 2, . . . , 30. The minimum of each curve is defined as the localization of the 
resonant state. 



In figure [5] we show the first thirty eigenvalues of the N = 100 p-block of the Hamiltonian 
matrix and the resonant energies £ calculated with the DO method. The curves D n (X) for 
n — 2, ... ,30 are shown in figure El Note the qualitative differences between the DO curves 
for I = and / = 1 in figures [3] and M respectively. These differences are due to the existence 
of a virtual state between the bound and the resonant states for the / = case, which it is 
absent in the I = 1 case, where the bound state is continued directly in a resonant state. 



In figure [7] we show the exact curve T vs. A and the approximate values obtained with 
two different basis-set size: iV = 100, for n — 2, . . . , 40 and N = 500, for n — 2, . . . , 200. As 
in the case with zero angular momentum, we obtain an excellent agreement between exact 
and approximate results. 
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FIG. 7. (color online) Exact (red line), N = 100 (turquoise dots) and N = 500 (black squares) 
values of T against the barrier height A for the I = 1 block of the potential Eq. (|34p . 



IV. CONCLUSIONS 

In this work we used a probabilistic interpretation of resonant (Siegert) states based on 
the conservation of the number of particles inside a time-dependent volume |8l-ll0|. The 
advantage of the probabilistic interpretation of resonant states is that is possible to work 
with Siegert states in a similar way to bound states, calculating probabilities, expectation 
values, etc. 

In particular, we obtain the exact equation (|22|) . which reduces to Eq. (125]) and Eq. fl37|) 
for / = and I = 1 respectively. These equations relate the inverse lifetime with other real 
magnitudes, the energy and the density of a resonance. In previous papers, once the energy 
of a resonance is obtained applying a real algebra stabilization-like method, the resonance 
width is calculated performing a fitting of the density of states ( 19| and references therein, 



201]) . In the present work Eq. (1221) gives a value for T with the same degree of accuracy 



that we obtain for the energy of the resonance £. 

We emphasize the simplicity of the calculations compared to other methods that use 



14 



complex algebra to study resonant states. We present our results for potentials with finite 
support, but Eq. ([22]) is valid in general, and the calculation of T could be corrected by a 
systematic perturbative expansion. 

An open question is if Eq. ( 1221) could be generalized for few-particle systems, where many 
channels are present. We are working in this direction. 
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